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ABSTRACT 


The  accuracy  of  proposed  lower  confidence  limits  for  system 

availability  is  analyzed.     Random  values  of  the  lower   100(1-  cO% 

A 

confidence  limit  Ac  _  .        .  for  system  availability  are  computed  for 

a  system  whose  failure  density  is  exponential  (  f\  )  and  whose  repair 
density  is  exponential  (/£.).      The   system  is  modeled  as  an  alternating 

renewal  process.      The   lOO^-o^    )th  percentile  point  of  the  generated 

A 

distribution  of  A„     .     /     ,  is   compared  with  system  availability  as  a 

A 

measure  of  accuracy  for  A„  ,  ,    ,    x. 

y  S  L(C<   ) 
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I.     INTRODUCTION 

The  system  considered  in  this  thesis  is  of  the  type  which  operate 
for  a  time,    fail,    are   repaired  and  returned  to  operation.      The   system 
then  may  be  said  to  have  two  states  and  can  be  modeled  as  a  renewal 
process. 

The  objectives  of  this  paper  were  the  following: 

1.  To  compute  exact  values  of  the  steady  state  availability  for 
a  system  with  given  distributions  of  times  to  failure  and  times  to 
repair. 

2.  To  compare  the  availability  of  a  system  based  on  the  test 
procedures  used  with  the  100(1-  &/  )  th  percentile  point  of  the  dis- 
tribution of  a  proposed  100(1  -  o(.  )%  lower  confidence  limit.  This 
will  allow  an  accuracy  analysis  of  a  proposed  lower  confidence  limit. 

The  method  used  was  to  fit  a  normal  distribution  to  simulated 
data  and  compute  the   100(1  -  o{  )  th  percentile  point  of  the  distribution. 

An  alternative  method  was  investigated,    which  fitted  a  2  para- 
meter gamma  distribution.     The  procedure  proved  to  be  inaccurate. 


II.       STATEMENT  OF  THE  PROBLEM  AND  ANALYSIS 

A.      DEFINITION  OF  AVAILABILITY 

In  this  paper  the  term  availability  [l]  is  defined  as  the  probability 
that  a  system  will  be  "up"  or  operable  when  called  upon  at  a  random 
time. 

Normally,    availability  is  not  specified  alone.      Both  reliability 
and  maintainability  are  usually  specified  together  with  availability. 
For  durable,    continuous  operated  hardware,    reliability  can  be  specified 
by  mean-time-to-failure  (MTTF).      Given  the  exponential  form  of  time 
to  failure  distribution,    which  is  applicable  to  most  types  of  complex 
equipment,    MTTF  is  a  constant,    the  reciprocal  of  the  failure  rate. 

Maintainability  can  be   specified  by  mean-time-to- repair  (MTTR) 
and  mean-preventive -maintenance -time.      Corrective  maintenance 
times  are  often  found  to  be  well  described  by  the  exponential  distribu- 
tion.    Given  that  the  time-to- repair  is  distributed  exponential,    MTTR 
is  a  constant,    the  reciprocal  of  the  repair  rate.     Accordingly  we 
make  the  following  definitions. 

(1)  A(t)     -    Availability  at  time  t,    or  the  probability  that 

an  item  will  be  operable  at  a  stated  instant  in 
time .     t  2  0 

(2)  A_  -        Interval  availability  is  the  time  average  of  A(t) 

during  intervals  of  length  T,    and  is  readily 
obtained  from  the  equation  for  the  average  value 
of  a  function. 


AT  — ^= -     ^      /      A(t)  dt 

/      dt  x    J 

c 

Alternately,    it  may  be  considered  that  there  is 

some  probability  distribution  h(t)  on  demand  time. 
Then  the  interval  availability  is  given  by 

DO 

E  [A(t)]  =     /    A(t)  h(t)  dt 

o 
In  the  special  case  where  h(t)  is  uniform  on  the 

interval  [0,T] 
h(t)    =    ^  0  £  t  £T 

o  t  y  t 

T 


\ 


thenAj,       =    ^,  J      A(t)  dt. 
o 
(3)    A    -  the   steady-state  availability  is  the  limiting  interval 

availability  as  T         ■>     o*z> 
T 


if 


A  =    lim    ■=,     /      A(t)  dt 
T->oO 


It  can  also  be  shown  that 


A     =    limA(t) 
T->   to 

If  the  components  do  not  undergo  checkout  or  repair  its  interval  avail- 
ability approaches  zero.      For  systems  with  exponentially  distributed 
failure  times  where 

A(t)=        e  lT 

^T  At 

A     ~>      A->    0  asT->    cO 


Thus  the  steady  state  or  long  term  uptime   ratio  of  a  component  which 
does  not  undergo  repair  approaches  zero. 

If  the  component  is   subject  to  repair,    its  steady  state  availability 
is  not  zero.     Assume  that  both  failure  and  repair  times  are  exponentially 
distributed  with  means    —    and    —      respectively.      Given  that  A(0)  =   1 
i.  e.  ,    component  available  at  start  of  mission,    then 

A,t)  =  -^—  +        X     .   c-t<A+  X  » 


JU^~     repair  rate 
A 


failure  rate 
then       A  =  lim  A(t)    =     — ■/■ -, — 

A^-t  A 

MTTF 


MTTF  +  MTTR 


Catron  [2]  has  shown  that  exact  values  of  A(t)  can  be  and  were 
computed  for  repair  distributions  other  than  exponential.     The  use 

of  interval  availability  as  t— *>  r>c    i.  e.  , 

T 

AT     =    I      J    A(t)  dt  -?»   A 

to  approximate  the  true  availability  over  the  first  mission  time  as  is 
currently  being  done  in  practice,    is  a  conservative  procedure. 

The  amount  by  which  the  limiting  interval  availability  under- 
estimates the  true  availability  during  the  first  mission  time  is 
dependent  upon  the  limiting  value,    i.  e.  ,    the  lower  limiting  value, 
the  more  it  underestimates  the  true  availability. 

(4)    A,-.     -    System  availability.      For  a  series  system  the 
system  availability  is  defined  as: 


K 

Ac    =    IT      A. 


where  A.  is  the  availability  of  the  component  of 
type  i.     i  =   1,  2,  •  • •  K 


B.      PROBLEM 

The  problem  considered  in  this  paper  is  to  check  the  accuracy  of 

A 

a  proposed  lower  confidence  limit  procedure    Ac  _  .       .  for  A     based 

upon  component  failure  time  data  and  repair  time  data. 

(1)  When  assuming  exponential  failure  rate     Cs.  . 

and  exponential  repair  rate     (O  .  then  A.  is   given  by 

A.     =       

The  following  ranges  of  the  failure  rate  (  C\  ■ )  and  repair  rate 
(  fO.)  has  been  used  in  the  simulation. 

0.  005  ±r<;  £    0.  01 
3.  Of   A^IO.  0 

which  implies  that  the  ranges  of  MTTF  of  component    i, 
denoted  by   A  .  and  MTTR  of  component  i,    denoted  by  j^L*. 
are  as  follows. 

100  £    A,  £  *z   200 

1 0  -  /^  -  3 

(2)  The  equation  of  component  availability  can  be  written  as 

a.   =    —&—  -—!—-  L_ 

1         /&+o(.L       i  +^       i  +  oicUu 

but    A.  =    -T—i- =   1   -0Q  tti  +  (oCLui,)Z  -  [tt£  lc-)3+  ■  •  • 


1 


TTK—  -  i  -y.  ♦  <y 


clM, 

8 


provided  the  series  converges,    which  it  does   since    oC^ JUL^\. 
Since  the  product   &£i  LLi   is  very  small,    in  this  case 
0.  0005±  *U  u.,,  6  0.  00333-  .  . 


A;  =     i  +w . —  =  l  -  etc  U 


K 


ThenAs  =  f  /?•   -,  7T  (l  -  U  U4*  l)  ==    \~    *2-<£lf**l 


Since  o<J  and  itk  will  have  small  variance  we  shall  apply  th( 
central  limit  theorem,    and  fit  a  normal  distribution  to 

k: 


S 
A  £ 

Z 


ThenAg    =     1  -   ^ott JJL± 


£  ln±ll    <*i  A* 


Assuming  that  A     ^  N  ]  A    ,      .^J  I     q.  & 


C.      ALTERNATIVE  METHOD 

An  alternative  method  was  investigated,    which  fitted  a  2  para- 


meter gamma  distribution  by  method  of  moments  to  the  distribution  of 

as-TT  k  -IT  jr—f- 
A        A  j 

whe re ^V    .  and  /?>  .  are  the  estimates  of   ,X  •  and    A,   .   respectively. 

Under  the  assumption  that  time  to  failure  and  time  to  repair  of 
the  components  were  distributed  exponentially  with  failure  rate  (X.   ■ 
and  repair  rate  lO    .   respectively.     The  following  approach  was  used 
to  make  the  fit: 


. Las-  i-i^(-^^-'i--^~ ~^) 


Q.  =  1-A.  =  1  -    . 
1  i 


OC^ 


die*  fie 

Thus  by  Taylor's  expansion 

-   ^Ac  =y_~  lU  (1-Q.)  =  T     (Q.  +    Q  i 
o      ' —  i         f —         i        — ^ — 


C-l 


2 


a.  and  b.  are  choosen  so  that 
l  l 

A 


^  /c:  a  k'  £       ,      /-)  x  \ 


E( 


X)  -  y, 


and 


A 

Q.  =  1 
i 


A. 


1  - 


A 


/A 


where  ^  .  and  /3   .  are  the  estimate  of  <X   .  and    /o  .   respectively. 
The  2  parameter  gamma,    I      {'C  ,\     )  was  fitted  by  method  of 
moments  to  the  distribution  of      j^  • 

E,    %  )  = 


J2- 
V 


r 


Therefore 


A 


E     ( 


y.)  _ 


Then 


Vcvr  (    y  ) 
A 

r  = 


r 


10 


s 


/v  /* 

by  replacing  all  oC  .  and    /3 .  by  CXi     •  and     /}  .. 
ThenZy^y.     is  distributed      ^'2 

_  C 

As  ■  e.  r 

From  this  assumption  a  one  sided  lower  100(1   -  o{  )%  C.  L.    for  A 
was  computed. 

A,  ^v* 
i.e..  P(2Y7iZJkt2r)  =  1-   <*£ 

p(as>  ^sl  c^c )  =  i  -  <*: 

r      ^    i  ^ 

where  [2  p"  J    =  smallest  integer  greater  than  or  equal  to  2  r 

The  procedure  showed  to  be  inaccurate  due  to  small  values  of    f"~ 


D.      PROCEDURE 

£   *     A. 

The  procedure  used  to  fit  a  normal  distribution  to  1  -    Z.  ^  i.  Za, 

as  described  in  B  is  based  on  computer  simulation. 

The  computations  of  the  estimates  of C\  •  i     and   1      -\     (1  =  1,  2,  *  *  "k] 
are  based  on  the  assumption  that  we  are  given  a  series  system  with 
K  types  of  components. 

V..     -    time  to  failure  of  the  j-th  component  of  type  i. 

The  time  to  failure  is  assumed  exponential  with 

failure  rate  C\  -• 

i 

T..     -    time  to  repair  of  the  j-th  component  of  type  i.     The  time 


to  repair  is  assumed  exponential  with  repa 


ir  rate  (6  .. 


V..  and  T..  are  independent. 
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(1)    A  random  sample  of  size  N.  components  of  type  i  are  tested 
until  failure  or  a  specified  planned  test  time  (PTT). 

T    . .  is  the  PTT  of  the  j-th    component  of  type  i 


a       if  the  j-th  component  of  type  i  fail,    i.  e.  , 
-*•       life  time  less  than  T    .. 


X.. 


j-      if  the  j-th  component  of  type  i  lives  beyond 

\J  HP 

oij' 


2__  X  i,  \       -    #  of  failures  of  type  i  components 


The  components  were  tested  against  the  (PTT)  T    ..  and  the  observed 

oij 

operating  times  V    ..  set  equal  to  V..  or  T    ..  whichever  is  the  smaller. 

oij  ij  oij 

(2)  A  random  sample  of  size  n.  of  failed  components  of  type  i 

are  repaired  with  repair  times  T. , ,    T.ot    ...,    T. 

r  r  il         lZ  ] 

(3)  Using  the  test  data  described  in  C  1  and  2  the  estimates  of 

Cv  .  and /L  were  computed. 
1      /      l  N- 

£.=   rxcj  hm 


in. 

i 


J*± 


JU-i 


Hi 


A 

(4)    With  the  data  in  C- 3  1000  values  of  A  i    .  were  computed. 

A      Ai      X 
where  A     =  1   -      2L    °^  C  L-L,C  and 

L-l      / 
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Z         is  the  100(1-  q(J)  th  percentile  point  of  the  N(0,  1)  distribution. 

E.      ANALYSIS 

A 
The  generated  values  of  A„  .  was  ordered  from  lowest  to 

highest  and  the   1000(  1  -  <X  )%  one  was  compared  to  Aq. 

A 

The  Ac      .       >  was  also  compared  to  A      and  the  probability  of 
o  i-/(oc. )  S 

success  was   computed  by  the  formula 

P  (ASL(<X     )   ^    AS) 
based  on  the   1000  values  of  A        .      .    ,. 

The   results  are  shown  in  Appendix  A  for  the  different  combinations 

of   oC    .,    /^  .,     N.,    n.  and  T    ...     To  avoid  fluctuations  due  to  the 
^    1     /«->  1  1         1  oil 

random  number  generator,    the  values  of  A_  •    .  at  different^    's 

o  -Li  (  u\.    ) 

was  computed  by  only  changing  Z^ 


13 


III.       ACCURACY  RESULTS 

By  changing  the  values  of  the  input  parameters  the  accuracy  of  the 
procedure  was   checked.     The  following  conclusions  were  made  based 
on  the  simulation: 

a.  Use  of  the   steady  state  availability  of  the  system  (A„)  to  ap- 

A 

proximate  the  estimated  availability  A_,  T  .     ,   .  over  the  mission  time 

is  a  conservative  procedure. 

b.  The  amount  by  which  the  steady  state  availability  differ  from 
the  estimated  100(1-  C\    )th  percentile  point  is  very  small. 

Example 

Case  3   -  Appendix  A. 
With  the  given  parameters  the  systems  steady  state  availability  is 
computed,    giving 

jc      A 

Ag  =    77^ =  .9925 

The  computed  95th  percentile  point  of  the  distribution  of  A„  _  ,~    «_* 

is  equal  to  .9931 

Thus,  |AS-AS(950)       =     |.  9925  -.9931        =  =  °°06 

Similarly  the  90th  and  80th  percentile  points  of  the  distribution  of 

A  A 

A„      .    ,  .  and  A     _  .    ?.   are  .  9930  and  .  9929  respectively.     By  increasing 

the  number  of  items  which  are  life  tested  from  50  to  100  and  keeping 
the  other  parameters  constant,    we  get  the  result  given  by  case  5. 
That  is, 

A       =     .  9925  (as  before) 

A 

ASL(0.  05)  *  -9928 

14 


This  implies  that 


A 

AS  "  AS(    950) 


. 9925  -  . 9921 


.  0003 


A  A 

In  this  case  A        .„„>  and  A^.    „„„,  =  .  9926  which  implies  that 


AS  "  AS  L(oC    ) 


=     .  0001 


for  d    -  0.  1  and  o<  =  0.  2 


A 


The  true  levels  of  confidence  associated  with  A  .  as  a  lower 

o  l_i(  tX.     ) 

100(1-  (X    )  confidence  limit  for  A     are   given  in  columns   10  through  12 

A 

for   CK  -  •  05,    .  10,    .  20  respectively.      That  is  for  Case  3,    A^  T  .    n_  . 

is  an  90.  8%  lower  confidence  limit  for  A„   rather  than  95%  lower 

A 

confidence  limit.      This  is  a  measure  of  the  inaccuracy  of  A„  T  .    „_. 

b  Li{.  vo ) 

as  a  lower  95%  confidence  limit  for  A„  for  the  given  parameter  values 

A  b 

in  Case  3.      Likewise  A        .     ,  _  .  is   really  a   84.  5%  lower  confidence 

A 

limit  and  A        .    ?_.  is   really  a  76.  9%  lower  confidence  limit. 

This  variation  between  the  true  level  of  confidence  and  the  proposed 

level  of  confidence  (e.  g.  ,    90.  8%  vs.    95%)  are  due  to  the  very  small 

A 

variance  of  A_  _  .     *     .  for  the  cases   considered.     It  is  felt  that  for 
t>  .L(  C\    ) 

these  cases  the  quantities 

A 

A  A 

S         S(l-  cx,  ) 
are  better  measures  for  the  accuracy  of  this  procedure. 
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APPENDIX  B 


LEV^L 


in 


''AIM 


HATE    =    71046 


LOO 


1000 


500 

6(  C 


roco 


\  T  i    • '  <;  T  n  N     \  1  f  M  !  5  )  ,  o  E  T  A  (  !  5  )  •  V  (  L  5  ,  1  0  0  )  ,  X  (  1  5  ,  1  0  C  )  ,  T  ( 
[  ISItYHATC 15) !n( 15)^1  (if  UAS1  (  1000  )  ,AS2  (  1000  ) 
,-lEAL     PTT(  15,  100/1500*2.0/ 
|«4N=1  (.00 
K  =  15 

Z  05=  1.645 

Z 10 =1.282 
Z2C= 0.342 

I  X=9  3  A£'Z3 

f  C PMAT(2F10.5,2I 10) 

Mj     (5»IO0)     »LPM  f  >  ,-3ETA(  !  )  ,M  I  )  , -N  I  {  I  > 
CC  4T  IMJE 

X 1=0 .0  

x  7 = <y .  o " 

X  3=0.0 

WR1XE(^»5JQ0J    

T  <  /  /  /  ,  3  3  X  , 

FCKMAT (30X, 3F10 
AS  =  1  .0 

7  C  C  0     I  =  1 

i   S=AS*6FTA( 

DC    8C00      '=1 

DC    2000     f=l 

'  =.N!LI.) 

K'OO    J=l 
x  (  I    J  )  - ' "'  •  n 

C  mi  L      i    \NPU<  I  X,  I  V  ,YFL  ) 
I  X  =  I  Y 


15,60) 
,  'vS3(  1 


05 /A  1/50 

000) 


as  ha  ror>' 

4) 


3 X  ,  'A SHAT1 C ' t  3X, «  ASHAT20 • , /// ) 


1  )/( 

,MN 

,K 

i 


,LE  A(  I  )  *-BETA(  I  )  ) 


=-1 
,J) 


?ooc 


■ 

3  3  ?  3 


V  I  I  ,  J 
1  F  (  V  ( 
[  F  (  V  (  I  ,  J  ) 

can  inue 
ccn  r  INUI 
p  :    3333    I 

K  M  =  N I  {  I  ) 


0/< 

LT 

GT 


LFA(]  ) *  -LT 
PIT  (  1  ,J  )  ) 
P  T  T  (  1  ,  J  )  ) 


J  t  Y  F  L  ) 
X  (  I V  J  ) 


I 


V(!,J)=PTT(l,J) 


K 


I  Y  ,  Y  F  L  ) 


r:n    ^v*   j  =  i  ,m 

C  ML     PANOU  (IX 

MI,  J)=-l  .0/6ETAU  l*ALOG(YFU 

CENT  IN Or 

V  A  =0.0 

AYrU=   )  .  0 

OG    4CC0     I  =1  ,< 

S  Ij I  v  T  =     . 

S U ^'XX--T.  0 

S  UM  V  =  C  .  0 

fcA-.tii  I  )  

5C  00    J=l  ,i   N 
SLWXX=S'U^XX  +  X(  I  ,  J) 
Slr"V  =  S'j;'V+-V  (  T  ,  J) 
SU.VT  =  SUMT+PTT<  T  ,  J) 
r.CM  IMJF 
SUMTT=0.0 
N  M  =  N I  {  I  ) 

DC    5  5  55    J  =  1  ,r..\!         

SLV,TT=S.JMTT+T(  I  ,  J) 

CC\T  INUE 

X  N = N  <  I  ) 

X  M  =  N I  (  I  ) 

AHAT  (  \   I  =  SUVXX/SUMV*_<  2 

YHAT(  I  )  =  nU'    TT/XM 

VA  =  VA+\flAT(  I  )*YFAT{  I  )**2/SUMT 

AYHA  T  =  AYHAT+AHA.J  (I  )  £-YHAT(  I  ) 

iOOO    CCNT INUE 

ASHAT=l .O-AYHAT 
AS  1  (  M  )=ASHAf-S  L)R  T  <  V  A.  )  *  7  05 
AS2( y)=ASHAT-SQKT< VA)*Z10 
ASM  V  >  =ASHA  T-S  DK1  {  VA  )^I?n 
1 F(AS.GF  .ASL  (M)  )     xl- XI  +  1  .0 
[F(AS.GC.AS2(M)  )    X2=X2+1.0 
I  f  (  AS.  Of.  .Ao3(M)  )     Xi  =  Xjil.O 


5  VO 


5^r- 


0*XN)/(2.0*XN+l.0_) 
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LEVEL        18  RTSTTN-  DATE    =    71046  09/41/50 

LfiOCO  .CCM.IN'UE       -      

DC    88C0    L=lf999 

JJ=L+1 

DO    8  800    J  =  J  J, 1000 


IF(  ASK  LJ  .  LT  .AS  I  (  J  )  )     G'l    TO    800 
TFNP1=AS1 ( L  ) 
AS1(L)=ASH  J) 
ASK  J)=TEMP1 

_8C.O     I  E(  AS2JLJ  .  LT.AS2(  J)  »_.£C_LQ_8S0 

TFMP2=AS2<L) 

AS2< L)=AS2( J ) 

AS2(J)=TEMP2 
R«0  f  F(  AS3(L)  .  t.T.AS3(J  )  )  GO  TO  8800 

TEMP  2=AS3(L) 

AS3( L) =AS3( J) 

AS3(J)=TFMP3 
_8800  CCNTINUE       .  __  _______     __   .__ 

WRITE (6,600)     ( ASK  I  )  , AS2( I  )  ,  AS3(  I  )  ,  1  =  1 ,  10  00) 

fcRITF(  6,200  ) 
?00    KRMATI///  ,3  2X  ,  '  Al.rA(  I  )  •  ,?Y  ,  •  RET  A{  I  )  •  ,7X,  »N<  t)  •  ,  5X,  »NN(  I  )  '  ,  ///) 

:    .<?<:•<.    i  =  1 ,  k 

I  T  E I 6 , 3 00 )     ALF A(  I) , 8ETA ( I ) , N ( I ) , N I ( I ) 
99-Q    CCNTINUF 
3 n:;    FORMAT  (30X,2FL0.5,2U0) 

X  f  N  =  v  \ 
ppn-3i  =  xi  />;  ,-i 

PR0H2=X2/X »N 

I  :      ,,.-)     \StP^Oai, PRCB2, PRG03.MN 
400    FORMAT (///, 10Xf • AS=»  ,F7.4,4X,'  PRUB1= ' , F7.4,4X, •PRCfl2=« ,F7.4,4X, 
1  i  pun  V3  =  «  ,F7.4,4X,  •  HN=«  ,14) 
hR I T  F ( 6 , 700 )     -.  S 1  (  9  50  )  ,  A S 2  <  900 )  ,  A S 3  (  8 00  J 
70  0    FiJKiMAT(///ilOAi_,ASll<z50.1-=J  if-7-.A_i-4.Xi  l.AS2{9.00_.=  "  •  T  7  .  4  ,4X,  •  AS 3  (  3  00)  = 
,F  l.< 
STCP 
END 


18 


APPENDIX  C 

A 

Derivation  of  Mean  and  Variance  of  A, 


*    A 


i 


fie 

/ 

but       ~r   -    Lu^     =  MTTR 


The  MTTR  should  always  be  less  than  one  mission  unit,    which  implies 
that 

oix  Ai  <  *<  i 


thus 


2       ,A,.  ,3 


Ai-  if(X.,  ,    =i-*i/u+(o^:)    -(^:)J  + 

provided  the   series  converges  which  it  does   since      0£^  /^-c    *--   J— 
Since  the  product    C\i.  Lulxs  very  small 

A.  =  -_- 5 *     1-OtiLui 

1     1+  o(c/^  / 

Thus 


As=    // 


^|-«<£i^  =     |-Z^/^c 


Since    (A.   .  and   ($>  .  will  have  small  variance  we  shall  apply  the  central 
limit  theorem,    and  fit  a  normal  distribution  to 

Ac  =  1  -  X  o(<  U-^ 


Thus  A  JC    A 

Ls 


2.  c<iM-i 


/\  ^=/  ^ 


when   /^^  is  the  estimated  MTTR  and  oC    ■  is  the  estimated  failure 
rate.     Based  on  the  assumption  that  time  to  failure  and  time  to 
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repair  is  distributed  exponentially  we  can  find  the  mean  and  variance 

A 

of  A„. 

k(ac)  =  i  -ZE&)E(A )  =  I -Tot  ip--  =  Hs 

s  L  /  *./ 

/S"  A- 

By  using  the  equations  given  by  OD  29304  [3],    forC^    .  and  Yew  {oC    .) 

we  can  find  the  variance  for  A~ 

q/    .         #  of  failures  2  N^ 

#  of  test  time  2  N    +  1 

i 
A 

E(  X   i  )  =  OC 

A  r^  / 

V^r-  (  (X  . )    -      __r^ 

i  ^zr 

A/.  i 

where  "/*     /     •  ~  sum  of  all  planned  test  times  on  component 

J  i  (in  mission  units). 


LCI    -     |  of  repair  ti 


A 

E( 


Thus  Vct,~  (As)  =         VO/r     (1  -  X  <*l  fit  ) 

]/e*r(e(.Li*-L)  (1) 

But  E(<^    .)  =7^(^)4^ 


+  CXC 


?tlJ 
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£(/.')  -  far  //.  )  +/d 


V"*  ( Hr) 


oLi  \!    z     tl:  +  l  \       J:  ,?z 


Therefore     V*T  ( UiJUC  )  =  (y£~.  +  <*■  C  )(/^  "7^   )  '""  ^/^ 


c 


2-  -> 


=  *i£* .  Hit'  +dLih~t.  ^ » o^i/L? 


ru\     n-      /      n^ 


j 

which  implies  that 


vgk* 


_L     °(i/*->i      f]L  +1 


(oLlUl) 


TUi       n. 


-j 
j 

Thus  by  ( 1 ) 

A  K- 


*  2   ^  M-L   .  (Ziti. 


A  A 


Since  <X  .  and    U-.  is  approximately  equal  to  X    .  and /-<-     respectively, 


and 

n.  +  1 
i 


n. 

l 


^iasn.     — =*>     oo 


we  can  write 

^                         -  .  *T    A        A            *         4 
E(AS)  =     /  —  2.    °^  CJU^L I* 

V^r   (A„)    =     2_     ^J±l 

s  ,  ^ 

^; 

4* 
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which  implies  that 

A  .  A. 

Ag  r*+    N(Ag,     Vcir  (As))  =  N(Ag,    V LLr    (Ag)) 


A         /^   Z 


-  n(i-  iS.,1^    z  &£} 


^  6SJ* 


Thus  by  the  above  assumptions  we  can  write 

7  (  ft*-/?**  £*W5)  =    '-  * 


c* 


SL('X  )      ~S  "  ^c< 


Nor  (/ts 


>c 


^        ^  2-      v     ± 


As'Z°<l*,    "75 


<*(L  y^o        \    Z. 
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